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Abstract. We consider the discretization of parabolic initial boundary value problems 
by finite element methods in space and a Runge-Kutta time stepping scheme. Order 
optimal a-priori error estimates are derived in an energy-norm under natural smoothness 
assumptions on the solution and without artihcial regularity conditions on the param¬ 
eters and the domain. The key steps in our analysis are the careful treatment of time 
derivatives in the H~^-noTm and the the use of an L^-projection in the error splitting 
instead of the Ritz projector. This allows us to restore the optimality of the estimates 
with respect to smoothness assumptions on the solution and to avoid artificial regularity 
requirements for the problem, usually needed for the analysis of the Ritz projector, which 
limit the applicability of previous work. The wider applicability of our results is illus¬ 
trated for two irregular problems, for which previous results can either not by applied or 
yield highly sub-optimal estimates. 
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1. Introduction 

Let H C be a bounded domain and T > 0 denote a time horizon. We consider the 
numerical solution of parabolic initial boundary value problems of the form 

u'{t) -|- A{t)u{t) = f{t), in 0 < t < T, (la) 

u{t) = 0, on dD, 0 < t < T, (lb) 

tt(0) = uq, in Q. (Ic) 

Here u' is the time derivative and Au = — div(aVu) -I- h ■ Vu -|- cu is a second order 
differential operator with coefficients a, 5, c G L°° and a{x,t) > a > 0, such that A{t) is 
uniformly elliptic for every point in time. It is well known, that for any uq G L^(H) and 
any / G L^(0,T; R“^(H)) there exists a unique weak solution in the energy-space 

W{t),T) = {u G : u' G L2(0, T; R-^(H))}, (2) 

which can be bounded a-priori by 

ll'“llw(0,T) = II^IIl2(0,T;R1(O)) + (3) 

< C'(||/||£,2(o,T;R'-l(n)) + ll^^o||L2(f2)). (4) 

Moreover, the constant C only depends on the domain and the bounds for the coefficients. 
We will refer to || • ||vi/(o,T) the energy-norm of the problem. If the coefficients and the 
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data are sufficiently smooth, the solution of (fTall - (fTc]) can be expected to be more regular: 
For instance, if one assumes that a, 6, c € then 

\\'^'\\L^io,T-,H^{n)) + (5) 

< C'(||/'||j^2(o,T;H-l(f1)) + 11/(0) “ ^(0)'WollL2(fI))> 

whenever the right hand side is bounded. If, in addition, also the domain is sufficiently 
regular and uq € Hq{^), then 

lkllL2(0,T;H2(Q)) + ||w'||L2(0,T;L2(n)) < C'(ll/llL2(o,r;L2(n)) + ||^01|//l(Q)) • (6) 

We refer to [9] for details and proofs of these estimates and further results. 

Let us emphasize at this point that the two estimates ([5|) and ([6]) also characterize the 
basic regularity spaces for parabolic problems. In particular, time derivatives typically 
lack two order of spatial regularity compared to the solution itself. 

In this paper, we study the numerical solution of (lla|) - (fTc]) by finite element discretiza¬ 
tion in space and Runge-Kutta time stepping schemes. Such approximations have been 
investigated intensively in the literature, see e.g. [g [13 H [ISl [531ES]; we refer to |5D] 
for a comprehensive collection of results and further references. Our main goal here is to 
derive order optimal a-priori error estimates in the energy-norm || • ||vf(o,t)- These can 
be obtained under minimal regularity assumptions on the coefficients and the domain, 
and we only require natural smoothness conditions for the solution. To keep the notation 
simple, we consider here in detail only the lowest order approximation by piecewise linear 
finite elements in space and the backward Euler method in time. The generalization of 
our arguments to approximations of higher order is however straight forward. 

To put our results into perspective, let us shortly recall some of the standard results 
for the Galerkin approximation of parabolic problems [2011531155] and compare them with 
the main contributions of our manuscript: 

As a first step, we will consider the semi-discretization in space by piecewise linear finite 
elements. The standard lowest order error estimate for parabolic problems then reads 

11^^ “ Uh\\L°°{o,T-,L^{n)) < Ch[\\uo\\Hi(n) + ll'w'liLqo.r (7) 

Here Uh is the finite element approximation and h the meshsize. This estimate amounts 
to [201 Theorem 1.2] with r = 1. Assuming sufficient regularity of the domain, the param¬ 
eters, and the solution, also second order convergence can be obtained here. The basic 
step in the proof of this error estimate is the decomposition of the error into 

u-Uh = {u-Rhu) + {RhU-Uh), (8) 

where Rh denotes the Ritz projection associated to the operator A; see [531 [25]. The proof 
of the error bound then relies on a discrete energy estimate for the semi-discrete problem 
and certain properties of the Ritz projection, in particular, on a bound 

||n - Rhu\\L2^n) < Ch\\u\\Hi{n), (9) 

for the L^-norm error, which can be obtained by the usual duality arguments mm- 
Without some additional regularity assumptions on the parameters and the domain, the 
estimate Q for the Ritz projection is however not valid, and therefore the validity of the 
bound ([7]) for irregular problems cannot be granted by the proofs given in [5U1 [531 [53] . Also 
note that the estimate ([7]) is somwhat sub-optimal concerning the regularity requirements 
for the solution: In fact, it should suffice to assume that u € L°°(0, T; which is 

already valid if u G L^(0, T; 71^(12)) n 71^(0, T; L^(n)). Let us emphasize that no spatial 
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regularity for the time derivative seems to be required. Also, in view of the a-priori 
estimate dG]), this latter condition would be a natural regularity assumption. 

As a replacement for ([7]), we will derive an error estimate in the energy-norm || • ||w(o,r)) 
which arises in the a-priori estimate for the solution. We will show that 

11"*^ “ Uh\\w('d,T) = 11^^ - 

< C'/i(||u||x,2(0,T;L/2(O)) + \W\\L^(0,T-,L^{n)))- (10) 

By continuous embedding, one has \\u —Uh\\L°°(o^T-,n) < C\\u — Uh\\w{o,T)^ which yields the 
corresponding estimate also for the error in the L^-norm at every point in time. One can 
easily see that the estimate m is optimal with respect to the approximation properties 
of the finite element spaces and also with respect to the smoothness requirements on the 
solution, which are natural for the problem under investigation. Also note that we require 
no additional regularity of the domain or the parameters for the proof of this result. 

As a second step, we will then investigate the time discretization by the backward Euler 
method. The standard error estimate for the fully discrete approximation reads 

Q™ax^||n(t”') -u^\\L‘^(n) < C/i(||no||//i(n) + \\ut\\L^(o,T-,H^{n)) +(H) 

see |20l Theorem 1.5] with r = 1. Here t = is the time-step size, and 

denotes the n-th iterate of the Euler method. Again, we observe a certain sub-optimality 
concerning the regularity requirments for the solution. Moreover, the proof given in m, 
see also [251 [23], is only valid under additional restrictive regularity assumptions on the 
parameters and the domain, which strongly limit the applicability of the results. 

As a replacement for m, we will derive the energy-norm estimate 

\W — Uh\\w{0,T) = 11^ “ ^/i||l 2(0,T;L/1(O)) + 11^^ “ '“h 11^2(o,T;L/-l(O)) (12) 

< C'(/l||u|||2(o,T;L/2(n)) + ^ll^1li2(0,T;L2(n)) 

+ T\W\\L^(o,T;H'^{n)) + t\W\\l^{o,t-,h-^ (n))] ■ 

Here Uh denotes the function obtained from by piecewise linear interpolation in time. 
The bound for again follows easily by embedding of W{0,T) into L°°(0, T; L^(H)). As 
before, the estimate (m can be seen to be optimal with respect to the approximation 
properties of the discretization, and, in view of the a-priori bounds ([5]) and ([6]), the regu¬ 
larity requirements for the solution seem natural. Moreover, no artificial regularity of the 
domain and the coefficients will be required for the proof of this estimate. 

One key step in the derivation of our results will be the careful estimation of time 
derivatives in the TT'^-norm. This seems natural in view of the definition of the energy- 
norm and its role in the a-priori estimates. The importance of such estimates has already 
been observed in the context of a-posteriori error estimation |15] . Let us also mention 
|131126] . where iL“^-Galerkin methods for the solution of parabolic initial boundary value 
problems have been investigated. 

A second difference to previous investigations is, that we use here a somewhat different 
error splitting as usually employed, namely 

U-Uh = {u-TThu) + {TrhU-Uh), (13) 

where is the L^-projection onto the finite element space. The L^-projection nh has 
important advantages in comparison to the Ritz-projector Rh- First, the estimate corre¬ 
sponding to do]) and further approximation properties can be proven without regularity 
assumptions on the domain or the parameters. In addition, tt/j commutes with the time 
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derivative, even for problems with time dependent parameters. At several places in our 
analysis, we will rely on the i?^-stability of the L^-projection mm- Morover, we use ap¬ 
proximation error estimates for the L^-projection in various norms, in particular including 
estimates in the ff“^-norm. 

The remainder of the manuscript is organized as follows: In Section [2l we formally 
present the problem to be investigated and we discuss the basic assumptions on the domain 
and the coefficients. In Section [3l we introduce the finite element spaces and summarize 
some estimates for the L^-projection that are required later in our analysis. Section [J] 
is then concerned with the derivation of the estimate m for the semi-discretization. In 
Section [5j we investigate the time discretization by the backward Euler method and we 
derive the second estimate (HI. For illustration of the wider applicability of our results, 
we discuss in Section [6] two irregular test proplems, for which, due to lack of regularity, the 
standard bounds © and m cannot be applied directly. Due to the weaker requirements 
for our estimates cni) and m, the optimal convergence in the energy-norm can however 
still be guaranteed also theoretically. We close with a short discussion of our results and 
highlight some possibilities for generalization and future investigations. 

2. Preliminaries 

Let D C M'^, d = 2,3, be some bounded domain and T > 0. We use standard notation 
for function spaces, see e.g. [9]; in particular, Hq{Q) denotes the sub-space of functions in 
H^{Q) with vanishing traces on dQ, and is the space of bounded linear 

functional on Hq{Q)-, the corresponding duality product is denoted by {-, •) 

Of particular importance for our analysis is the energy-space 

W{0,T) = {n G : u' € L\0,T-, (14) 

which is equipped with the norm ||u||iy(o,T) = II^^IIl2(o,T;H1(o)) + lk'llL2(o,T;H-i(n))- For 
later reference, let us recall the following embedding result 0 Sec. 5.9]. 

Lemma 1. Let u € W{0,T). Then u ^ C{[0,T]; and 

sup ||u(t)||i2(f^) < C'||u||iy(0,T)- 
0<t<T 

As a consequence, all estimates derived in the energy-norm || • ||u/(o,T) automatically 
yield corresponding bounds in || • \ i-®-; pointwise in time. Here and below, we 

denote by C some generic constant which may have different values in different occasions. 

Let us now turn to the initial-boundary value problem (fTaD - (fTc|l . We assume that the 
operator A has the form Au = — div(aVti) -|- b ■ Vu -|- cu. In order to proof well-posedness 
of the initial boundary value problem, we assume that the parameters satisfy 

(Al) o,c € X (0,T)) and b € x {0,T)f, and 

(A2) a{x,t) > a for some constant a > 0 and a.e. {x,t) € D x (0,T). 

Because of the first assumption, A defines a bounded linear operator from L^(0, T ; Hq (D)) 
to For a.e. t G (0,T) we then define an associated bilinear form by 

a{u,v;t) = / a{x,t)\/u{x)\/v{x) + b{x,t) ■ Vu{x)v{x) + c{x,t)u{x)v{x) dx. (15) 

Jn 

The weak formulation of the initial boundary value problem (|lap - (fTc] ) now reads 
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Problem 2 (Weak formulation). 

Given f G L'^{0,T] and uq G L‘^{Vt), find u G W(0,T) such that u(0) = uq and 

(n'(t),u)j^-i(f^)x_H'i(n) + a{u{t),v;t) = (/(t), r’)//-i(n)x_H'i(r 2 ) 
holds for all test functions v G Hq{GL) and for a.e. t G (0,r). 

Under assumptions (A1)-(A2), Problem [2] has a unique solution u G 1U(0,T) and there 
holds ||u||vt/(o,T) < C'(II'Wo||l 2 (o) + ll/llL 2 (o,T;H-i{n)) with a constant C that depends only 
on the bounds for the coefficients and on the domain; see [9] for details. The proof relies 
on a Gronwall argument and the following properties of the bilinear form. 

Lemma 3 (Continuity and Carding inequality). 

Let (Al) and (A2) hold. Then there exist constants Ca,a,r] > 0 such that 

a(u,'(;;t) < C'a||u||j^i(n)||'y||//i(n) and a{u,u;t) + r]\\u\\l 2 ^^^ > a\\u\\j^^^^, (16) 

and these estimates hold uniformly for a.e. t G (0,T) and all n, u G H^{Ll). 

Proof. The continuity follows from the Cauchy-Schwarz inequality. Using the lower and 
upper bounds for the coefficients, we get 

a{u,u-,t) > a||Vu||| 2 (f^) + ||&(i)||L-||Vu||i 2 (f^) ||u||i 2 (o) + \\c{t)\\L--\H\l 2 (^Qy 

The Carding inequality then follows by Young’s inequality and choosing the coefficients, 
for instance, as a = a/2 and rj = a/2 + || 6 ||^oo/( 2 a) + ||c||loo. □ 

3. Properties of the L^-projection onto finite element spaces 

For the semi-discretization in space, we will employ a standard finite element method. 
To avoid technical difficulties, we assume that Ll is polyhedral and that it can be partitioned 
into a set Th = {T} of simplicial elements T. More precisely, we require that 
(A3) 7h is a regular simplicial partition of fl, i.e., the intersection of two different ele¬ 
ments either empty, or a vertex, an entire edge, (an entire face) of both elements; 
(A4) Th is locally quasi-uniform, i.e., there exists a 7 > 0 such that the diameter hx of 
an element T and the radius px of the largest ball that can be inscribed in T are 
related by 'yhx < Pt ^ hx for all elements T. 

Civen such a mesh Th, we consider the standard finite element space 

Vh = {v€ : v\x G Pi(r) for all T G Th}, (17) 

of piecewise linear continuous functions that vanish on the boundary. Furthermore, we 
denote by iTh ■ L‘^{Ll) —>■ I 4 the L^-orthogonal projection defined by 

{TThU,Vh)L2{n) = {u,Vh)L2{n) for all G 14- 

Obviously, || 7 r/iu| 42 (Q) < ||ri||L 2 (Q), i.e., the L^-projection is stable (a bounded linear op¬ 
erator) on Under the assumption (A3)-(A4), it is however also stable on 

Lemma 4 (77^-stability). Let (A3)-(A/) hold. Then 

hhu\\m{n) < C\\u\\Hi{n) for all u G 

and the constant C only depends on the domain and the regularity constants of the mesh. 

For globally quasi-uniform meshes, the result follows directly from the Bramble-Hilbert 
Lemma and an inverse inequality. The proof for locally quasi-uniform meshes has been 
given in [4]; see [3] for generalizations including higher order approximations. 

We will also require the following approximation error estimates. 
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Lemma 5. Let (A3)-(A4j hold. Then 

11 ''^ - -^hnWH^in) < Ch’^~^\\u\\jjk^^^ 
for all u G Hq with 0 < k <2 and —1 < s < min{l, k}. 

Proof. For completeness, we sketch the main steps. The case s = 0 and 0 < A; < 2 is well 
known and follows from the Bramble-Hilbert lemma and scaling arguments. To show the 
estimate for s = 1 and 1 < /c < 2, let us denote by : Hq{LI) —>■ I 4 the ^-orthogonal 
projection defined by 

{7r’fu,Vh)HAn) = {u,Vh)Hpn) for all Vh G 14. 

Recall that \\tt\u — u 11/11(0) < C'h’^ ^11 ^11//* (o) for t 4 k < 2, which is the usual finite 
element error estimate [5]. We can then proceed by 

\\u - 7rhu\\Hi{n) < 11?^ - T^hTrinWH^n) + hh{u - TTiu)\\Hi(n) 

< (1 + C)\\u — < (1 + C)C'h\\u\\H2(Q), 

where we used the projection property and the //^-stability of tt/i, and the approximation 
properties of vr^ in the last two steps. Now assume that u G with 0 < A: < 2. Then 

- vrfeu||//-i(o) = sup (u - 7r/iU,u)i2(n)/||u||//i(o) 

(O) 

= sup {u,v-7rhv)L2(^n)/\HHAn) <Ch’^"^^\\u\\Hk^Q). 
veH^in) 

Here we used the approximation property of tt/j for s = 0. This yields the estimate for 
s = — 1 and 0 < k < 2 and completes the proof. □ 

Corresponding estimates for real valued s and k can be obtained by interpolation argu¬ 
ments. Also note that the dependence on h in the above estimates can be localized. 

4. Galerkin semi-discretization 

As a first step in the approximation process, let us investigate the semi-discretization 
in space by finite elements. Proceeding in a standard fashion, we define 

Problem 6 (Semi-discretization). 

Find Ufi G H^{0,T-,Vh) such that Ufi{0) = T^hUo o,nd 

{u'h{t),Vh)L^Q)+a{uh{t),VK]t) = (/(A),fh)//-i(o)x//i{n) (18) 

for all test functions Vh G 14 and a.e. AG (0,T). 

The first term could also be written as iih)//-i(o)x//i(n)' By choosing a basis for 

Vh, the semi-discrete problem yields an ordinary differential equation, and existence and 
uniqueness follow by the Picard-Lindelof theorem. More precicely, we have 

Lemma 7 (Well-posedness and discrete a-priori estimate). 

Let (Al)-(A4) hold. Then Prohlem\^has a unique solution Uh G fA^(0,T;14) and 

\\Uh\\wiO,T) < C{\\uo\\l2^q'^ + ||/||l2(0,T;//-1(O)))- 

The constant C in this estimate depends only on the bounds for the coefficients, the domain, 
and the constants characterizing the regularity of the mesh. 
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Proof. It remains to verify the a-priori estimate. By testing the discrete variational prob¬ 
lem in the usual way with Vh = Uh{t), one can show that 

< vlWWh^n) + ^II/WIIh-i(o)- 

Here we only used the Carding inequality (11611 and some elementary manipulations. Inte¬ 
grating over time and applying a Gronwall argument then yields the energy estimate 

ll“lli°°(o,r;L2(o)) + \\'^\\h(o,T-,H^{n)) < 11/11^2(0 ,T-,H-pn)))- 

To obtain the remaining bound for the time derivative, observe that the fl~^-norm of a 
discrete function can be expressed as 

\\u'h{t)\\H-pn) = sup K(t),'c)L2(n)/||u||Hl(0) 
v£H^{n) 

= sup «(t),7r/ju)i2(Q)/||'(;||j^i(n). 

v£H^{Q) 

Using the discrete variational problem with Vh = T^hV further yields 

< {\\fhit)\\H-'^{n) + C\\uhit)\\Hi{n))\\'^hv\\Hpn)- 

Since tt/i is bounded on Hq{Q), we further have \\'n'hv\\H''^{n) < C'||u||j:^i(f 2 ). Using these 
estimates in the expression for the norm of the time derivative finally yields 

\Wh{t)\\H-'^{n) < C{\\fh{t)\\H-i(n) + lkfe(^)llL7i(f^))' 

The result then follows by integration over time and the energy estimate derived before. □ 

As a direct consequence of the definition and the similarity of the continuous and the 
semi-discrete variational problems, we obtain 

Lemma 8 (Galerkin orthogonality). 

Let u and Uh denote the solutions of Problems andl^ respectively. Then 

{u{t) - Uhit),Vh)H-i{n)xH^^{n) +aiu{t) - Uh{t),Vh]t) = 0 

for almost every t € (0,T) and all Vh G 14. Moreover, 7r/j(n(0) — Uh{0)) = 0. 

We can now turn to the error analysis of the semi-discretization. To this end, we divide 
the error into an approximation error and a discrete error as 

u{t) - Uh{t) = {u{t) - TThUit)) + {'KhU{t) - Uh{t)) = (i) + (a). 

Due to the approximation properties of tt/i, the first term can be bounded readily by 

Lemma 9 (Approximation error). Let (A3)-(Af) hold. Then 

11^ - 7r/iU||w(0,T) < C'^(lkllL2(0,r;L72(O) + lk^llL2(o,T;L2(o))), 
and C only depends on the domain and the shape-regularity constants of the mesh. 

Using Galerkin orthogonality and the discrete stability of the method, the discrete error 
component can then be bounded as usual by the approximation error as well. 

Lemma 10 (Discrete error). Let (A1)-(A4) hold. Then 

WlThU - Uh\\w{ 0 ,T) < C\\u — 'n'hu\\w{ 0 ,T)- 

Moreover, the constant C only depends on the domain, the bounds for the coefficients, and 
the shape-regularity of the mesh. 
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Proof. The discrete error eh = T^hu{t) — Uh{t) satisfies 

{&h{t),Vh)L^(n) + a{eh{t),Vh]t) = {TThu'{t) - u'{t),Vh)L^(n) + a{T^hu{t) - u{t),Vh]t) 

= a{Trhu{t) — u{t),Vh; t) =: {f {t),Vh) 

and also e/i(0) = T^huit) — Uh{t) = 0. Using the continuity of the bilinear form, we have 

ll/(^)llLI-l(n) < Ca\\u{t) - 'Khu{t)\\Hl(n)- 

By application of Lemma [7] to the equation for the discrete error, we thus obtain 
l|eft||ty(0,T) < C'll/llL2(0,r;_H'-l(f1)) - - Vr/jU||vi/(0,T); 

and this already proves the assertion of the lemma. □ 

By combining the previous two lemmas, we readily obtain our first main result. 

Theorem 11 (Energy-norm error estimate for the semi-discretization). 

Let (Al)-(A4) hold, and letu anduh denote the solutions of Problems\^and\^ respectively. 
Then 

Ik - 'U.h\\w{0,T) < C'^(lkllL2(o,r;fl2(r2) + 11^11^2(0,T;L2(0))), 
and the constant C only depends on the domain, the bounds for the parameters, and the 
shape regularity of the mesh. 

Let us emphasize that the regularity requirements for the solution are natural. In view 
of approximation properties of finite elements, they are also almost necessary to guarantee 
the approximation order one in the energy-norm. Also note that no additional regularity 
of the domain or the coefficients was required for our proofs. 

5. Time stepping 

Let us now turn to the time discretization of the semi-discrete Problem [6] by the back¬ 
ward Euler method. For ease of presentation, we only consider here uniform time steps of 
size T = T /N, and therefore set = nr. To allow for evaluation of the coefficients and 
the data at individual points in time, we further assume that 
(A5) a,6,c € iyi’°°(0,r;L°°(Ll)), and 
(A6) f 

The condition (A5) could be relaxed by using a Discontinuous-Galerkin method for the 
time discretization. In addition, we require that the time step is sufficiently small, i.e., 
(A7) T < I/t], 

where r/ is the constant from the Carding inequality (|15l) . This condition could be avoided 
by treating the lower order terms in the bilinear form in an explicit manner. To facilitate 
the notation, we will use in the sequel 

to denote the discrete time derivatives at Applying the backward Euler scheme for 

the time discretization of Problem El then leads to 

Problem 12 (Fully discrete scheme). 

Set u\ := TT/iUo and find G 14 for 1 < n < N , such that 

{dru'^+^,Vh)L^n)+a{ul+^,Vh;t^^^) = {f{t'"^^),Vh)H-pn)xH^{n) 
holds for all test functions Vh & Vh- 


( 19 ) 
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The equation (fT9l) amounts to an implicit time-stepping scheme and that the snapshots 
can be computed recursively. The assumption (A7) guarantees that the elliptic problem 
for each time-step is uniquely solvable. More precisely, we have 

Lemma 13 (Discrete well-posedness and a-priori estimates). 

Let (Al)-(A7) hold. Then Problem,[IE has a unique solution {u'j^}o<n<N, o,nd 

- C'(||'Uo||l 2 (q) -h 

The norm on the left hand side is a discrete version of the energy-norm || • ||w(o,T)- 
Proof. Testing (fT9|) with and proceeding as in the proof of Lemma [7| yields 

<^ii<‘Siiy(n)+»ii“r'ii!-(n)+^ii/(<’‘*')iiL.(n)' 

Via a Gronwall argument, we then obtain the discrete energy estimate 
maxi<„<jv \\uh\\h(n) 

Due to condition (A5), the last term can be estimated via Talyor expansion by 

— ^(ll/lli2(o,T;L/-l(r2)) +''"^11/ llL2(o,r;LI-l(r2))) • 

In order to derive the estimate for dru'f^, recall that 

Wdrul'^^Wn-iin) = sup 7r,ju)i,2(Q)/||u||j^i(n). 

Using Vh = T^hV as a test function in the discrete scheme (fT^ . we further obtain 

= (/(*”+'), 

< C(||/(C+i)||„-.(n| +CK+'||«.(n,)lkft.^|l„.(n)- 

Using the stability of the projection tt/j, this allows to estimate the discrete time derivative 
by known terms. The assertion then follows from the previous estimates. □ 


Let us now turn to the derivation of error estimates for the full discretization defined 
in Problem m Similar as in the previous section, we use an error decomposition into an 
approximation error and a discrete error by 

u{n -ul = {u(f) - + {'KhUif ) - ul) = (i) + {in). 

The hrst component can be bounded by the approximation estimates for vr/j as follows. 

Lemma 14 (Approximation error). Let (A3)-(Af) hold. Then 

'^n=l “ drTThU{n\\H-i^^-^ + ||n(r) - TThU{t'^)\\H^n)) 

— ^(^^ll^lli 2 ( 0 ,T;Ll 2 (Q)) + h‘^\\u |li 2 ( 0 ,T;L 2 (n)) llL 2 (o,r;LIl{ 0 )))' 

Proof. Using Taylor expansion in time and the properties of the projection, we obtain 
'^^T\\drU{r) - dr7rhUit^)\\Hi^^^ 

— ^11^ “ '^hU lli2(0,T;H-l(n)) — ^ h^\\u lli2(o,r;L2(r2))' 
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In a similar manner, we obtain for the second term 

< Cdlu — Tfeu||^2(0,T;//l(n)) + ''”^11'“ “ '^hU llL2(o,T;//l(n))) 

< C"(/l^||u||2,2(o,T;H2(fl)) + T‘^\W\\L2(0,T-,H^{n))) ■ 

The assertion then follows by summing up the two contributions. □ 

The discrete stability of the scheme (|19l) allows us to bound the discrete error as follows. 
Lemma 15 (Discrete error). Let (A1)-(A7) hold. Then 

Yl^n=l 

— ^(^^ll^llL2(0,r;iL2(n)) +''■^11^ llL2(0,T;LI'l(n)) llL2(0,r;LI-l(L2))) ' 

Proof. Let e)) = drT^hu{t'^) — denote the discrete error. Then 
'9r(e)(+d Uh)i2(Q) + a(e”+d u;,; r+^) 

= {dru{r^^) - u'(r+^),Ufe)i2(Q) + a(7r,,u(r+^) - u(r+^), u,*; r+^). 

With similar arguments as in the proof of Lemma [T3l we then obtain 
max IIII r 2 (-O') 

l<n<A ^ 


I n||2 

\^h II ^2/'0^ \ / 


’■lleSllHi(n) 


< CY.^T(\\dMn - !.'(«”)+ \\n.n(e) - 
Using Taylor expansion, the hrst term on the right hand side can be further bounded by 
''^^T\\drU{f') - < ^||w^dlL2(0,T;LI-l(n))’ 

In a similar way, we obtain for the second term 

< C^W'KhU — 'w|Il2(o^T;H 1(^)) ^ llL2(0,r;Hl(^))) 

< C''(h^||li|li2(4n^in + l.J^2(f^)) + \W\\\2 (^■a ^^■a + l . 

Here we employed Taylor expansions and the approximation and stability properties of 
the projection tt^. By combination of the two estimates we obtain 


max^ Ikhu(r) - <||^2(t2) + ^^T||7rfeu(t”) - 


— ^(^^ll^lli2(o,r;_H-2(n)) +"^^11^ lli2(o,T;L/i(^)) 


+ T n‘ 


"\\2 


IIl2(o, 


T;LI'-l(n)) 


)• 


Using the characterization of the H ^-norm for finite element functions, the dehnition of 
the discrete scheme dUD, Galerkin-orthogonality, and some basic manipulations, we can 
further bound the discrete time derivative terms by 


V T\\drTThU{t'^) 

^ 


- Uu 


l|2 


+ ||„(t") - w(t")|||,.,n| + IK - ir,KC)||„.(n))' 

The first and second term on the right hand side can be estimated as above, and the third 
term is an approximation error which has already been bounded. □ 
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Summing up the two estimates for approximation error and the discrete error yields 
Lemma 16 (Discrete energy-norm error estimate). Let (Al)-(A7) hold. Then 

< C'(/l^||u|||2(0,T;H2(n)) + ^^II“^IIl2(o,T;L2(o)) 

++''‘^11^^ lli2{0,T;_f/-l(n))) ■ 

To finally obtain an estimate in the energy-norm || • ||vk(o,t)) l^t us denote by 

t-t^ 


^n+l _ ^ 

^h{t) = - -u'f, + 


^n+1 _' ^»^+l _ ’ 

the linear interpolant of the fully discrete approximations in time. Using the dehnition 
of Uh and simple manipulations, we now obtain the second main result of this paper. 


n+l 


r <t<t 


n+l 


( 20 ) 


Theorem 17 (Error estimate for the full discretization). 

Let (Al)-(A7) hold and u be sufficiently smooth. Then 

\\u - Uh\\w(0,T) < C'(/l^||M||i 2 (o,T;H 2 (n)) + \\u'\\L^(0,T-,L‘^{n)) 

+ "^^11^ lli2(0,T;Lll(f^)) +''"^11^ 11+2(0,r;LI-l(L2)))' 

Like in the previous estimates, the constant C here only depends on the domain, on the 
bounds for the coefficients, and on the shape regularity of the mesh. 

Let us emphasize that in view of the bounds (l5|) and Q, the regularity requirements 
for the solution are reasonable and that the rates are optimal with respect to the ap¬ 
proximation properties of the discretization. Also note, that no additional smoothness 
assumptions for the domain or on the spatial regularity of the parameters were required. 


6. Numerical tests 

For illustration of the benehts of our results, let us shortly present two irregular test 
problems, for which the standard error estimates © and m cannot be applied directly, 
while our estiamtes (|10l) and (1121) still provide order optimal error estimates. 

6.1. Lack of smoothness. Let us consider the one-dimensional heat equation 

dtu{x, t) = Uxx{x, t) -|- /(x, f), 0 < X < 1, 0 < t < 1, 

with homogeneous initial and boundary conditions. Here dtu = u' is used for the time 
derivative synonymously. We assume that the exact solution has the form 

u(x,t) = y^ Un sin(n7rx) sin(n^7rt) with = (1-|- 
^ —^n>l 

The parameter e > 0 is assumed to be small. The special form of the solution allows to 
compute various norms of u analytically and, in particular, to show that 

ue L‘^{0,T-,H‘^{Q)) and u'£ L‘^{0,T; L‘^{Q)). 

At the same time, one can easily check that 

L°°{0,T-,H^+^^n)) and u'^ L\0,T-, 

Due to lack of regularity, the standard a-priori estimate d?]) for the semi-discretization 
therefore cannot be applied, in contrast to our energy-norm estimate (|1U|) . which yields 

||u - tt/i||L°o(0,T;L2(o,7r)) < C'll^ “ '>J‘h\\w{ 0 ,T) < C'h. 
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Note that, since u ^ -L°°(0, T; (Q,)), the rate of convergence for the error in the norm 

of L°°{0,T; L'^{Q)) cannot be improved substantially here. For a solution of the form 

u(x,t) = 'S^ Un sinfnvrx) sin(n^/^7rt) with = (1 + 

— ^n>l 

one can see in a similar manner that 

ue L‘^{0,T;H^{Q)), u'e L‘^{0,T; and u" G ^^(O, T; 

However, u" ^ L^(0, T; L^(n)), and therefore the standard estimate (fTTI) for the full dis¬ 
cretization cannot be applied due to lack of regularity. On the other hand, our error 
estimate m for the full discretization still allows to guarantee 

||u(r) - u/j(r)||i 2 (!^) < ||u - Uh||Loo(i 2 (Q)) < Cllu - Uh\\w(0,T) < C{h + r). 

Note that by interpolation one has at least u G L°°(0, T; here, so the rate of conver¬ 

gence for the error ||tt(t) — Uh{t)\\i 2 (^Q-^ in terms of the meshsize may possibly be improved. 
For sufficiently smooth solutions, the estimates of [231 El EH would in fact predict the 
optimal rate \\u{t) — u/i(t)||i 2 (Q) < C{h‘^ + r) here. 

6.2. Discontinuous parameters. As a second test case, we consider a thermal diffusion 
problem on a square covered by an inhomogeneous medium. The governing system reads 

dtu{x,t) = div{a{x)'Vu{x,t)) + f{x,t), xG(— 1,1)^, 0 < t < 1/2. 

As before, we presecribe homogeneous initial and boundary conditions. Moreover, we 
assume that the diffusion parameter has the form 

I 1, Xi ■ X2 > 0, 
a[x) = < 

xi ■ X2 < 0 , 

where e is some small positive constant. It is well-known [IS]> that the associated elliptic 
operator Lu = — div(aVtt) for such a problem is rather irregular. More precisely: for 
every /3 > 0, one can choose an e > 0, such that L is not an isomorphism from 
to (n). In particular, the maximal value of /3, such that 

\\RhU - u\\L2(n) < Ch^\\u\\fji(^^) ( 21 ) 

holds for arbitrary u G can be made arbitrarily small. Note that the standard esti¬ 

mate ([7]) does not apply directly here, since only (I7TD holds instead of Q. A generalization 
of the error estimate to non-smooth problems however still allows to guarantee 

11^ “'“/i||L°°(0,T;L2(n)) < (22) 

provided that the solution u is sufficiently smooth; see 1201 Sec 19] for details. Since /3 can 
be aribtrarily small in general, this estimate is highly unsatisfactory. Our energy-norm 
estimate (uni) for the semi-discrete approximation however still applies and yields 

||t6 - U/i||l°°( 0 ,T;L 2 (q)) < C\\u - Uh\\w{0,T) < C'h, 
provided that the solution has the required smoothness. This results hold regardless of 
the spatial regularity of the diffusion parameter o('). For illustration of the validity of this 
theoretical result, we also provide some results of numerical tests. To verify the conver¬ 
gence rate for the semi-discretization, we integrated the semi-discrete problem numerically 
in time with a very accurate time stepping scheme. The resulting errors obtained on a 
sequence of uniformly refined meshes are summarized in Table [TJ As predicted by the 
theory, we observe convergence of the energy-norm error with hrst order. The numerical 
tests actually yield a better convergence rate for the error in the norm of L°°(0, T ; L‘^{Q)), 
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h 

0.50000 

0.25000 

0.12500 

0.06250 

0.03125 

rate 

ei 

1.65010 

0.82062 

0.39475 

0.19318 

0.09710 

1.03 

62 

0.37632 

0.11211 

0.02928 

0.00741 

0.00186 

1.93 


Table 1. Errors ei = ||ii - Uh\\w{o,T) and 62 = \\u - Uh\\L^(o^T-,L^in)) 
obtained on a sequence of uniformly refined meshes with meshsize h. 


which in fact is the optimal one from an approximation point of view. We however cannot 
give a full explanation for this observation yet. 

7. Summary 

Various results concerning the numerical analysis of Galerkin approximations for par¬ 
abolic problems are available in the literature; see [20] for a comprehensive overview and 
further references. This paper contributes to this active field with providing error esti¬ 
mates in the energy-norm ||u||vi/(o,T) = ll'i^llL2(o,T;Hi) + seems to be a 

natural choice from an analytic point of view, but which has not been studied intensively 
in previous works. 

In this manuscript, we considered only low order discretizations of a simple model 
problem. The general approach is however applicable to much wider class of problems and 
discretization schemes. Morover, our arguments may also be fruiteful for the derivation 
of a-priori error estimates in other norms da [13 [23 [Ml ES] and for the derivation of 
a-posteriori error estimates [3 El 113 [ 13 . 

Our main motivation to consider the energy-norm, was to overcome a sub-optimality 
of the standard estimates m and m concerning the regularity requirements for the 
solution, for the domain, and for the parameters. This sub-optimality is partly due to a 
loose handling of time derivatives in the estimates, and, on the other hand, stems from the 
use of the Ritz projector in the error decomposition, which requires duality arguments and 
regularity of the underlying elliptic problem. We therefore utilize here the L^-projection 
in our error splitting and carefully estimate time derivatives in the 77~^-norm. 

In our presentation, we focused on a-priori error estimates in the energy-norm || • ||vu(o,T)) 
and we could establish optimal convergence rates under minimal regularity assumptions. 
By continuous embedding, we could also obtain estimates pointwise in time with the same 
convergence rates. In our numerical experiments, we observed a better convergence of 
the error in the norm of L°°(0,T;L^) for a particular problem, which does not follow 
directly from our results. For sufficiently regular problems, this better convergence is 
well explained by the standard results |23[23[23]. A justification for the case of certain 
irregular problems is however missing. 
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